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lO ■ Abstract: In the case of large 1-3 mixing angle as sin 2^13 > 0.03, we investigate the 

possibility for measuring the leptonic CP phase by using only v^ -^ u^ oscillations inde- 
pendently of v^ -^ fe oscillations. As the result, we find that the best setup to measure 

t^^ , the CP phase without strongly depending on the uncertainties of the other parameters is 

^-^ ' around the energy E = 0.43GeV and the baseline length L = 5000km. In this region, 

the CP phase effect remains even after averaging over the neutrino energy. We also find 

^ ' that there is a CP sensitivity even in the short baseline length L < 1000km if Am-gj^ is 

C^ • determined with an uncertainty of about 0.1%. In the T2KK experiment, we explore the 

possibility for measuring the CP phase by using a baseline from Tokai to Kamioka after 
determining A?7i|^ at the baseline to Korea. As the result, we find that some information 
of the CP phase can be obtained from the measurements. 
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1. Introduction 

The finite mass of neutrinos and the mixings among different flavors have been confirmed 
in various neutrino experiments. Recently, the LSND anomaly [1] is refused by the result 
of MiniBOONE experiment [2] and the results obtained in the past neutrino experiments 
can be almost explained by the neutrino oscillations among three generations. The effect 
of such neutrino oscillations has been considered not only in the experiments on the earth 
but also in the far universe for example supernova explosions [3, 4, 5]. 

For the values of mass squared differences and mixing angles, the results of atmospheric 
neutrino experiments [6], K2K experiment [7] and MINOS experiment [8] provide 

lAmiil ~ 2.5 X 10"^eV^ sin^ 2(923 ~ 1 (1-1) 

and the results of solar neutrino experiments [9] and KamLAND experiment [10] provide 

Am^i ~ 8.1 X 10"^eV^ sin^ 6^12 ~ 0.31. (1.2) 

On the other hand, only the upper bound for the 1-3 mixing angle 

sin2 26li3 < 0.16 (1.3) 

is obtained from the CHOOZ experiment [11]. We cannot determine the sign of Amg^ at 
present from the experimental data. Furthermore, we have no information on the leptonic 
CP phase 5. Unveiling these unknown parameters is one of the most important aims in 
the next generation neutrino experiments. In particular, the value of 5 is very important 
at the view point of the leptogenesis [12]. 



One of the serious obstacles in determining the value of 6 is the eight-fold degeneracies 
[13, 14, 15]. There are some cases in which the uncertainty of 5 becomes large due to 
the effect of degeneracy. One of the turning points is whether ^13 can be determined 
by the next generation reactor experiments like Double CHOOZ experiment [16] and the 
superbeam experiments like T2K experiment [17], NOz/A experiment [18]. In this paper, 
we concentrate on the cases that sin 2^13 is larger than 0.03 and assume that ^13 will be 
found in the Double CHOOZ experiment without being affected by the 6-5 ambiguity. In 
this case, it is suggested by many authors that the remaining four-fold degeneracies can be 
also solved within a few decades. One possibility is the observation of neutrinos from the 
same beam to two detectors on different baselines like Tokai-to-Kamioka-Korea (T2KK) 
proposal [19] and SuperNOz/A proposal [20]. The other possibilities are the combination of 
more than two kinds of neutrino sources, long baseline plus atmospheric neutrinos [21, 22] 
and long baseline plus reactor neutrinos [23]. The possibility for using the Wide Band 
Beam and analyzing the spectral information of neutrino events is also investigated in ref. 
[24]. In these proposals, the measurement of the leptonic CP phase is also explored by 
using v^ -^ Ve (i^e -^ T^fj) osciUatious. 

In our previous work [25], we suggested that we can measure the leptonic CP phase 
by using only v^ -^ v^ oscillations in the region E < 2GeV, L > 2000km if sin^ 2^13 
is large. In the analysis, the probability for i/^ -^ f^ oscillations can change about 0.4 
by the CP phase effect and then there remains a difference of about 0.2 between the 
maximal and minimal values of the probabilities even after averaging over the neutrino 
energy. In these considerations, we have explored new possibilities of experiments to be 
performed after a decade in addition to solving parameter degeneracies and determining 
the value of the CP phase. It is easy to observe the v^ -^ u^ oscillations in superbeam 
experiments and neutrino factory experiments. There is a potential to measure the CP 
phase 5 independently of v^ -^ v^ oscillations in the case that the parameters except for 
5 can be measured precisely. If we assume the unitarity in three generations, the channel 
of u^ -^ Vf. oscillations should be related to that of z/^ -^ v^ oscillations. This means 
that we can predict the behavior of one oscillation channel by the behavior of the other 
channel. (See references [26] about the discussion of the unitarity in the lepton sector, for 
example.) If we find the difference between this prediction and the experimental result, we 
must consider effect of new physics and we will have some constraints to the unified theory 
in the high energy physics. Thus, the measurements of b in two independent channels are 
very important for exploring the new physics beyond the Standard model. Recently, the 
exact formulation of neutrino oscillation probabilities has been extended to the case of non- 
standard interaction in view of the era of precision measurement of parameters [27] . A lot of 
investigations have been also performed about the exploration of non-standard interaction 
and the non-unitary effect by future experiments. See [28], [29] and the references therein. 

In this paper, we give the detailed analysis for the possibility of the measurement of 
the CP phase in v^ -^ z/^ oscillations, considering the systematic error, uncertainties of 
parameters except for b and the uncertainty of matter density, which were not considered 
in our previous work [25] . We investigate the best conditions of baseline length and energy 
for measuring the CP phase in v^ -^ z/^ oscillations by using both analytical and numerical 



methods. We also calculate the baseline dependence and ^13 dependence of CP sensitivity. 
Furthermore, we consider the reason why we have no CP sensitivity in comparatively short 
baseline length and show how the CP sensitivity depends on the uncertainty of Amg^. 
We estimate how this uncertainty can be improved in future experiments and how CP 
sensitivity can become good in comparatively short baseline experiments. We usually use 
the v^ -^ v^ channel in order to measure Am|j^ precisely. Namely, we need two different 
baselines to measure the CP phase in the case of short baseline by using this channel. So, 
we consider the T2KK experiment as concrete setup. We use one baseline from Tokai to 
Korea to measure the value of Amg^ as precise as possible and use another baseline from 
Tokai to Kamioka to measure the CP phase. 

Outline of this paper is the following. In section 2, we explore the energy and the 
baseline length where the CP dependence in v^ -^ v^ oscillations becomes large. In section 

3, we assume the concrete experimental setup and calculate the CP sensitivity. In section 

4, we clarify what is the problem for measuring the CP phase in relatively short baseline. 
Then, in section 5, we reinvestigate the CP sensitivity in the case that this problem can 
be improved. In section 6, we conclude. Finally, in appendix, the approximate formula for 
the coefficient of cos 5 is derived, which is applicable for the case of large L/E. 

2. Region with Large CP Dependence 

In this section, let us review the CP dependence of neutrino oscillation probabilities ob- 
served in the superbeam experiments. We also investigate the energy and baseline length 
where the CP dependence becomes large. 

The Hamiltonian in matter is represented as 

H = 02^TH'T^Ol^, (2.1) 

where O23 is the rotation matrix between the second and the third generations and F is 
the phase matrix defined by F = diag(l, l,e ). Without loss of generality, we can factor 
out the part of ^23 and 5. This makes it possible to include the matter effect only in the 
reduced Hamiltonian H' as 

H' = Oi3Oi2diag(0, A21, A31 jOfaOfg + diag(a, 0,0). (2.2) 

Due to this factorization, we can transparently understand how oscillation probabilities de- 
pend on the CP phase 5. In the above expression, we use the equalities Ajj = Amfj/{2E) = 
{mf - m'j)/{2E), a = V2GFNe ~ 7.56 x 10"^ • pYe, Gp is the Fermi constant, Ne is the 
electron number density, p is the matter density, Yf, is the fraction of electrons, E is the 
neutrino energy and rrii is the mass of ui. The CP dependences of the probabilities for 
i^iJL -^ ^e oscillation, v^ -^ v^ oscillation and t'^ -^ v-r oscillation are given by 

Pfjie = Afj^e COS 5 + B^e s'lU 6 + Cfj,e, (2.3) 

Pf^f,^ Af,f, COS 6 +C^^, (2.4) 

P^r ^ B^rsm6 + Cf,r, (2.5) 



in ref. [30]. Eqs. (2.4) and (2.5) hold exactly in the case of ^23 = 45° [31], where the 
coefficients ^^q, .B^q and C^q (ck = e, /i, r) are the quantities determined by the parameters 
except for 6. If we assume the unitarity in the framework of three generations, the sum of 
these probabilities has to be one for any value of S. This leads to the relations among the 
coefficients 

V + ^M^ - 0' (2-6) 

Bf,e + B^r = 0, (2.7) 

C^e + C^^ + Cf^r = 1- (2-8) 

It is well known that the order of magnitude for the coefficients is represented as 

^M/^ = -\e = Oiasis), B^r = - V = ^(asia), (2.9) 

C^e = 0(«2) + 0(s?3), C^^ = 0(1), C^r = 0(l), (2.10) 

by using the small parameters a = A21/A31 and sin^is = S13. See ref. [31] for example. It 
should be noted that the magnitude of A^^ is as same as that of A^^, and is not so small. 
Therefore, v^ -^ v^ oscillation will be one of the important channels in order to obtain 
the information on the CP phase, although it is hard to observe the CP phase compared 
to v^ -^ Ve channel because of the large CP independent term C^^. If we observe some 
differences between the values of the CP phase measured by the two independent channels, 
this means the violation of the unitarity in three generations and we can obtain some of 
the important clues for new physics. 

Below, it is considered how we should choose the energy and the baseline in order 
to measure the CP phase by using only u^ -^ v^ oscillations. At first, we numerically 
calculate the region in E-L plane with large |^^^|, which is the coefficient of cos 5. In 
this calculation, we use the following parameters, Am.3]^ = 2.5 x lO^^eV , sin 2^23 = Ij 
Am^i = 8.1 X 10-5eV^ sin^ 6^12 = 0.31, sin^ 26*13 = 0.16, p = 3.3g/cm3 and Y^ = 0.494. 

In figure 1, the black color shows the region with large |^^^|, namely at relatively low 
energy E < 2GeV and a long baseline L > 2000km. In other words, \A^^\ becomes large in 
the region with approximately L(km)/ii^(GeV) > 2000. However, in present experiments, 
only the event rate averaging over the energy can be observed in the region with such large 
L/E due to the finite energy resolution of the detector. Hence, we need to learn whether 
the CP phase effect remains in such situation and we would like to investigate the condition 
under which the CP phase effect becomes the largest after averaging. 

In order to investigate the behavior of vl^^ more accurately, we derived the approximate 
formula of A^^ from the exact one [32] in the Appendix. In the derivation, we kept in mind 
that our approximation should be valid for the energy region E ~ 0.1-lOGeV realized 
in superbeam experiments. More concretely, we defined Aj as the effective mass of i-th 
neutrino divided by {2E) and we took the approximation Ai < A2 ^ A3, a <C A3 and 
5^3 <^ 1. Then, we have left only leading order terms of small quantities, A21, Ai, A2 and 
S13. In order to be a good approximation for the region with large L/E, we did not neglect 
the term with the order of 0{A2i) included in the oscillating part and multiplied by L/E. 



IDGDD 



BDDD 



^□DD 



4DDD 



■SDDD 




;>0.07E 



Figure 1: Region with large lA^^^j. In the black region, the magnitude of \Af^^\ becomes large. 



Hence, the approximate formula for A^^ is calculated as 
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where J^ = S12C12S23C23S13C13, \j = Aj - Aj, A • ■ = AijL/2 and 
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In eq.(2.11), A^^ is represented as the sum of two terms Ai and ^2- ^1 is slowly oscillating 
term according to the change of energy as controlled by A21. A2 is rapidly oscillating 
term as controlled by A'^i- In the small L/E region, Ai can be neglected and the main 
contribution comes from ^2- As the value of L/E increases, Ai also gives the contribution 
and A2 oscillates faster. Therefore, only Ai remains in the region with sufficiently large 
L/E and after averaging over the energy. The total behavior of ^^^ can be described as 
the oscillation around the average value determined by Ai . The coefficient of sine function 
in Ai is given by 



4JrA2i(a-A2i cos 26*12) 

A2 
^21 



4JrAmli{2aE - Am^i cos 26*12) 
{2aE - Am\ cos 26*12)^ + Am\^ sin^ 26*12 ' 



(2.15) 



If we use the parameters sin^ 2^23 = 1 and sin^ ^12 = 0.31, it is found that the value of 
local maximum is given by 

,max _ 2 J^ ^.^2 / V2Am|^sin2^\ ^ sin26'i3 ^.^2 / V2Am|isin2^\ 



^ sin 26^12 \ 4Ei I 4 \ 4E, 

at the energy 

^ ^ Amji (cos 2^12+ sin 2^12) ^ ^ ^^^^^ _ Amj^ _ 3.3g/cm=^ 

^ 2a ' 8.1 X 10-5eV2 p ' ^ ' ' 

If we fix the energy at this value and from the maximal condition yzATTigi sin2^i2i/4£'£ = 
(2n + l)7r/2 (n = 0, 1, 2, • • • ) for (2.16), we can determine the baseline length Li as 



V2Ei(2n + l)TT , , , 8.1 X lO-^eV^ 

V4^ ^77^ = 5000km • 2n + 1 —^ 

Am^^sm 26^12 Am^^ S.Sg/cm"^ 



^ v--n--^-;" ^ 5000k^ . (2n + 1) . "-^"^^ ^^. (2.18) 



If we use the average density calculated in the PREM [33] corresponding to each baseline, 
Ai becomes maximal at Lg = 5000km and 10000km in the earth mantle. We also find from 
(2.18) that yl5"^^ attains to about 0.1 in the case of sin^ 26*13 = 0.16. 

Next, let us consider ^2- The factor (4JrA2i/A2i) sin A2X included in A2 does not 
change so much compared with sin(2A3]^ — A2x) for the change of energy. Namely, we 
regard 

4JrA21 . r, 2J,.Am^i . 

— = smA2i ~ ^sm(aL) (2.19) 

A21 o,E 

as the amplitude of A2 and roughly speaking, the magnitude decreases inversely propor- 
tional to the energy. Next, let us consider the oscillating term related to Ag^^. The construc- 
tive interference of Ai and A2 occurs in the case of sin(2A3]^ — Ag^) — sin(2A3i — a)' ~ —1. 
This leads the maximal condition {2Am'^-^ — 2aE)L/{AE) = {An — 1)tt/2. In other words, 
A^^ becomes maximal near the energy 
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En{GeV) = dv^km (n = l,2,3,---). (2.20) 

^ ' (4n-l)7r + 2-2-1.27ai|^ ^ ' ' ' ^ V ; 



From this expression, it is expected that the peak appears near E = 6.7, 2.9, 1.8 • • • GeV 
in the baseline L = 5000km. In the case of small energy, we can approximate as 

4 • 1 27^"^^^ — 
i?n(GeV) c. ^-^fpH^ (n = l,2,3,---), (2.21) 

[An — l)TT 

then the difference between the neighboring peaks is given by 





= En- 


1.27^i^ 3.0- 


Ami L 


AE(GeV) = 


2.5-10-3cV2 5000km 

n? 



(2.22) 

at the region satisfying the condition n ^ 1. Roughly speaking, we can observe the 
oscillation in the case that the energy resolution cte is smaller than the difference between 



the neighboring peaks. For example, the energy resolution is given by dg = 0.085GeV in the 
case of Water Cherenkov (WC) detector used in the Super-Kamiokande (SK) experiment. 
From the above condition AE' > (jg, we obtain the condition n < 6. This means that we can 
distinguish up to the fifth peak (E ~ IGeV) at the baseline length of L = 5000km. In the 
energy lower than E ~ IGeV, we cannot distinguish the different peaks and have to take 
the average. Namely, the effect of A2 can be neglected. In figure 2, left and right figures 
show the magnitude of A^^ and P^^ as the function of energy at the baseline L = 295km 
and 5000km . 

In the left figures, the blue and the red lines represent the magnitudes of ^^^ and Ai. 
In the right figures, the red and the blue lines are corresponding to the probabilities in the 
case of (5 = 0° and 5 = 180° respectively. One can see that the magnitude of A^^ is small 
in the case of L = 295km from the top left figure. On the other hand, the CP phase effect 
is large in the case of L = 5000km even after the averaging and the value becomes maximal 
around E = 0.43GeV. This coincides with the result obtained by the analytical expression. 
We also find that the peaks appear in the position calculated by A2 as we discussed before. 
In the middle right figure, we have large CP dependence also in the survival probability. 

Next, let us consider the case of inverted hierarchy. We obtain the coefficient of cos 6 in 
the case of inverted hierarchy by the replacement A31 -^ — A31 in (2.11) (exactly speaking, 
we perform the replacement A31 -^ — A31 + 2A21 cos 2^12). Ai does not change in this 
replacement. Therefore, the total event rates in both hierarchies are almost the same 
in the low energy region. On the other hand, A^^ of inverted hierarchy becomes slightly 
smaller in high energy region because of the MSW effect [34] due to the 1-3 mixing. Totally, 
the CP sensitivity is expected to be small in the case of inverted hierarchy compared with 
the case of normal hierarchy. 

Here, we comment the case of anti-neutrino oscillations. We obtain the probability 
Pfifi for the anti- neutrino oscillations i?^ -^ v^ by the replacement 5 — > — (5 and a ^ — a in 
P^^- As there exist only CP even terms in -P^/i, we only have to replace a to —a. We can 
obtain the energy E^, where Ai is maximal, as 

^ Ar.i,(- cos 2^,2 + sin2^,2) ^ ^^^^^^ Arn|^ _ 3^3gA^ 

k 8.1 X lO-^eV^ p ' ^ ' 

by the same procedure performed in P^^- However, Ei is nearly the production energy 
of /U and so it would be hard to observe the CP phase effect for this energy. So, we have 
to measure the CP phase in higher energy region than Ei for the case of anti-neutrinos. 
Considering also the smallness of the cross section of anti-neutrino, which is about 1/3 of 
that of neutrino, the use of neutrinos has an advantage compared with the use of anti- 
neutrinos for the measurement of the CP phase through i/^ (i/^) events. In figure 3, we 
show the energy dependence of Ap_p_ and P^^ in the baseline of L = 5000km in left and right 
figures. In the left figure, the blue and the red lines represent the value of Ajip_ and Ai. In 
the right figure, the red and the blue lines correspond to the probabilities in case of (5 = 0° 
and 5 = 180° as same as in figure 2. We can see that the absolute value of Ai attains up 
to 0.1 around the energy E = 0.1-0.2GeV. We also find that the sign of A^^ is negative 
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Figure 2: Energy dependence of A^^ and P^^- Left and right figures show the magnitude of A^^ 
and P^fj,. In the left figure, blue and red lines represent the magnitudes of A^^ and Ai respectively. 
In the right figure, red and blue lines correspond to the true values of (5 = 0° and 180°. 

over the entire region. This comes from the denominator of Ai, which has negative sign 
for the replacement of a — > —a 

Let us summarize the results obtained in this section. 

• In i/^ ^ v^ oscillations, the averaged value of A^^ becomes maximal around the 
energy E = 0.43GeV and the baseline length L = 5000km and 10000km in the earth 
mantle. 



The average value of Aau is positive in the region E > 0.2GeV. This means that we 
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Figure 3: Energy dependence of A^p and Pp,fi. Left and right figures show the magnitude of Ap_ji 
and Pfip,. Blue and red lines represent the values of A^^ and Ai respectively in left figure. On the 
other hand, red and blue lines in right figure represent the probabilities with 6 = 0° and 180° 

only have to measure the total rate of f^ events in order to observe the CP phase 
effect and we need not the good energy resolution of the detector. 

• We obtain more information on the CP phase by observing the energy dependence 
of f^ events at E > IGeV in addition to the total rate. 

• In relatively short baseline like L = 295km, the magnitude of A^^ becomes small and 
the observation of the CP phase effect is difficult unless other parameters except for 
S are precisely known. 

• In the case of inverted hierarchy, the CP phase effect included in the total rate is 
similar to the case of normal hierarchy. The CP dependence in the high energy region 
is slightly reduced. 

• Considering the production energy of fi and the cross section, it is easy to observe 
the CP phase in the experiment by neutrinos compared with anti-neutrinos. 

• As the CP phase effect is proportional to S13, it is difficult to observe in the case of 
small S13. 

In the above discussion, we do not consider the decrease of the flux of neutrinos accord- 
ing to the distance. From the statistical point of view, the short baseline is advantageous 
because of the large event numbers. However, in the case that the ratio of A^^ to the 
probability is small, the observation of 6 is strongly affected by the uncertainties of other 
parameters except for 6 as we will discuss later. 

3. Estimation of Signal from Leptonic CP Phase 

In this section, we estimate how precise the CP phase can be measured in v^ — > u^ os- 
cillations by using x^ method. As an experimental setup, we consider the 4MW and 2.5° 
Off- Axis JPARC beam and WC detector with the fiducial mass of 500kt [17], which is the 



same as those of T2HK experiment. We take two kinds of baseline length L = 295km 
and 5000km as in the previous section. We assume ten years data acquisition by using 
only neutrinos. We use the same parameters as in figure 1. We also assume that the mass 
hierarchy has already been determined before measuring of the CP phase in v^ -^ u^ oscil- 
lations. For example, the possibility of determining the mass hierarchy by the atmospheric 
neutrinos was discussed in ref. [21]. They concluded that the observations using the 545kt 
WC detector in three years determines the mass hierarchy at 2-a C.L. if sin^ 2^13 > 0.05. 
The parameter uncertainties, except for (5, are assumed as 5% for ^12, ^23 and ^13, and as 
4% for Am^i, and 1% for Aml;^ by expecting an improvement for the next ten years [35]. 
We also consider 5% uncertainty of matter density. The energy window for our analysis 
\s E = 0.4-1.2GeV and is divided into 20 bins. We calculate Ax^ by using the energy 
dependence of QE f^ events and the total rate of CC u^ events. As the backgrounds, we 
consider NC Ux events, where x takes all flavors. We set the uncertainties of signal and the 
background normalization rig, n^ as a{n'^^) = 2.5%, a{ns ) = cx), a{nh) = 20% and the 
uncertainties of their tilts tg, th as a{ts) = 2%, o-(tf,) = 2%. As the value of a{n^'^) = 2.5% 
may be optimistic, we discuss how the change of signal normalization affects the results 
later. The reason for setting a{ns ) = cx) is to be normalization free for QE events and to 
prevent the double counting in QE and CC events. We use common normalizations and 
tilts except for this. The expected number of QE events Nqei observed in the i-th bin and 
CC events Ncc are calculated as 



NQE^ = sf"^ (1 + nf^ + tr 
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(3.2) 
where s and b are the signal and background for the case that the uncertainties of normaliza- 
tion and tilt are not considered. See [36] as the more detail of the definitions. Furthermore, 
we define Ax^ as 
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by using NQEi and Ncc, where a in Ax^^^ represents the sum for QE and CC events and 
i in Axpara represents the sum for the parameters. Namely, Xi corresponds to the mixing 
angles, the mass squared differences and the matter density. As the true value of the CP 
phase, we take Stme = 90° and 180° and plot the value of Ax^ as the function of test value 
in figure 4. We use the Globes software [37] in this calculation. 
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calculations including only systematics and the blue including also parameter uncertainties. Left 
and right figures correspond to 6true = 90° and 180°. Top and bottom figures are calculated under 
the assumption of normal hierarchy in the baseline length L = 295km and 5000km. 



In figure 4, left and right figures correspond to the case of Stme = 90° and 180°. Top 
and bottom figures show the Ax^ calculated in the baseline L = 295kni and SOOOkni for 
the case of normal hierarchy. Prom the top figures, we can find that the measurement of 
the CP phase is difficult in L = 295km because of the uncertainties of parameters. On 
the other hand, the difficulties are largely decreased for the case of L = 5000km. We can 
find from figure 4 that the allowed range is 50° (106°) in l-a (2-0") C.L. for 6true = 90° 
and 88° (120°) in 1-a {2-a) C.L. for 6true = 180°. We do not show the figures for the 
case of inverted hierarchy, because they are almost similar to the case of normal hierarchy 
and only the CP sensitivity becomes slightly worse. This can be understood by using the 
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Figure 5: L dependence of CP sensitivity. In left and riglit figures, we take 5true = 90° and 180°. 
Red, green and blue lines show 1,2 and 3-ct C.L. lines respectively. We use the same parameters as 
in figure 4. 



approximate formula of A^^ given in (2.11) as follows. Namely, in low energy region, A^^ 
does not depend on the sign of Amg^^ because of the averaging of A2 and in high energy 
region, the number of events is suppressed due to the MSW effect compared to the case of 
normal hierarchy through A2. 

Next, we investigate how the CP sensitivity depends on the baseline length and the 
magnitude of ^13. The value of CP sensitivity is shown in figure 5 and figure 6. The 
CP sensitivity stands for how wide CP angles within 360° is allowed at certain C.L. when 
we set a true value of the CP phase. In the case that the allowed range is narrow, we 
can determine the CP phase precisely. Left and right figures correspond to the case of 
^true = 90° and 180°. Here, only the case for normal hierarchy is shown as we obtain 
similar result for inverted hierarchy. 

Figure 5 shows the baseline dependence of the CP sensitivity. In the both cases, the 
best sensitivity is realized around L = 4500 ^ 5000km. This coincides well with the result 
obtained in (2.18). In the baseline length of L = 10000km, the CP sensitivity is not good 
because of the small statistics due to the too long distance. 

Figure 6 shows the ^13 dependence of the CP sensitivity at L = 5000km. We can see 
that the CP sensitivity becomes worse gradually according to the decrease of 6*13 in both 
cases. This can be understood by eq.(2.11) as ^^^ is proportional to S13. It is also found 
that the CP sensitivity for 5true = 180° is good compared with 5true = 90° in 3-0" C.L. This 



is interpreted as follows. In f^ 



oscillations, the probability depends on the CP phase 



through cos 5. So, we obtain the largest differences between the probabilities for the two 
extremes, (5 = 0° and 5 = 180°, and the probability for 5 = 90° is in between these two 
cases. This reduces the difference to other probabilities and the allowed range becomes 
wide. 

Next, let us consider the dependence of CP sensitivity on the systematic errors. As the 
backgrounds in u^ disapearance channel are small compared to the signal, o"(nfe) and cj{ti,) 
hardly change the results. The main factor affecting the results is the signal normalization 
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Figure 6: 0i3 dependence of CP sensitivity a.t L — 5000km. In left and right figures, we take 
Strue = 90° and 180°. Red, green and blue lines show 1,2 and 3-a C.L. lines respectively. 
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Figure 7: a{n'^'-') dependence of CP sensitivity at i = 5000km. In left and right figures, we take 
Strue = 90° and 180°. Red, green and blue lines show 1,2 and 3-ct C.L. lines respectively. 



Figure 7 shows the dependence of the CP sensitivity on Uf^ beam normalization at 
L = SOOOkm. We can see that the CP sensitivity becomes worse gradually according to the 
increase of a{n'^ ) in both cases. In this baseline, the value of CP phase is mainly deter- 
mined by the total events averaged over the energy. As the total events are proportional 
to the cos 5 through Jj- in eq.(2.11), it is considered to be difficult to distinguish the effect 
of 6 and the signal normalization. 

In the end of this section, we consider the reason why the CP sensitivity is so bad for 
the case of L = 295km. For the case of the short baseline, we could in principle determine 
the value of 6 with high accuracy if the parameters, except 6, are precisely determined 
because of the large statistics. As the flux is inversely proportional to L^, the statistics 
for L = 295km becomes about 250 times larger than those for L = 5000km. However, 
the advantage is lost if the other parameters are including a few percent uncertainties. 
This is due to the small contribution of the CP phase effect to the probability and the CP 
phase effect is hidden by the large contribution from other parameters. The next question 
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Figure 8: (5(Ato|]^) dependence of CP sensitivity at L = 295kni. In left and riglit figures, we take 
^true. = 90° and 180°. Red, green and blue lines show 1,2 and 3-cr C.L. lines respectively. 

is which parameter prevents the determination of b. It is found that the uncertainty of 
Am|^ is related to the precise determination of b by the numerical calculation as suggested 
in [38]. In figure 8, we show the dependence of the CP sensitivity at L = 295km as the 
function of the uncertainty of Ayrig^ in the case of normal hierarchy. The experimental 
setup is assumed as the same as that in figure 4. The left and right figures correspond to 
the cases for btme = 90° and 180°. Horizontal axis represents the 5{Am'^i) / Am'^^ in the 
unit of %, and the vertical axis represents the value of CP sensitivity. From these figures, 
it is found that the CP sensitivity becomes gradually good below an uncertainty of 0.6%. 



4. Correlation between 6 and the Uncertainty of ^ml^ 

In the previous section, it was shown that the uncertainty of Atti^^ becomes the serious 
obstacle in determining 6 precisely. In this section, we investigate these correlation in more 
detail by using the analytical expression and explore the possibility for the improvement. 
For the case of 0(A2i) <C 1, the probability P^^ is reduced to the well known expression 



P ~ 



OTA 

1 — sin^ Ag;^ 1 sin A21 sin A'^^ cos Ag;^ cos 6, 



A 



(4.1) 



21 



as shown in appendix. If we perform the replacements A31 — > A3i(l + e) and 5 -^ 5' , the 
probability changes as 



P^^^l-sin^A;,! 



8J,.A2i 



1 - sin^ A' 



31 



A21 

8J,A2i 

A 



sin A21 cos 6' + A^ieL j sin Ag^^ cos Ag;^ 

A3ieA2iL 



sin A91 sin A01 cos Aoi cos 6' + 



^21 



^31 



31 



21 



8JrA2i sin Ag;^ 



(4.2) 
(4.3) 



up to the leading order of S13 and e. Here, we should note that the energy dependence 
of the correction term from the uncertainty of ATTig^ is approximately equal to that of 
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Figure 9: Correlation of 6 with the uncertainty of Amgj^. In left and right figures, we take 
-E = 10 GeV and 0.7 GeV. Regions with same color have almost the same probability given by 



baseline lengths 



PZl'^-Pp-i^l < 0-001 (0.00005) in the left (right) figure. Different colors correspond to the different 



the cos J' term. Namely, we cannot distinguish two probabilities P^^(A3i(l + e),6') and 
P^^{A3i,6) when the relation 



cos 



S' 



COS( 



A3ieA2i-L 
8Jr.A2i sin Aj 



~ cos S 



21 



Am'^^eaL 
SJrAml^sin' 



(4.4) 



is satisfied, where we use the approximation A21 — a. This means that the value of the 
CP phase cannot be determined if e is larger than a certain value. This is the correlation 
between 6 and the uncertainty of Am^i and is a serious obstacle for measuring the CP phase 
in v^ — > v^ oscillations. The small value of the CP sensitivity in the baseline L = 295km 
is due to this correlation. 

Next, let us estimate the magnitude of e giving the small CP sensitivity. If we substitute 
a = 7.56 • 10-^ pYe, p = S.Sg/cm^, Y^ = 0.494, J, = si2Ci2S23C23Si3cf3 ~ 0.045, Am|i = 
2.5 • lO-^eV^ and Am^^ = 8.1 • lO'^eV^ into (4.4), we obtain 

cos(5 ~cos(5-— — -— -— J— . 4.5) 

sm (3.1 • 10 *L) 

In the case of relatively short baseline, the above relation is further reduced to 

cos 6' ~ cos 6 — 170e (4.6) 

and does not depend on the baseline length. Let us illustrate the meaning of this relation 
in figure 9. 

In figure 9, regions with same color have almost the same probability as that for 
Strue = 90°. More concretely, the region with |P^^ - P^™^| < 0.001 (0.00005) is plotted in 
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the left (right) figure. The horizontal axis and vertical axis are taken as cos 5' and A?7i|^ 
(e) respectively. Left and right figures represent the case for E = lOGeV and O.TGeV. 
Red, blue and black colors correspond to the cases for L = 295km, 1000km and 5000km. 
It is obvious that the value of 6 can be determined by measuring the probabilities for two 
different energies for the case of L = 5000km from figure 9, because a superposition of 
the black curves for the two different energies (left and right figures) would lead to a clear 
intersection point representing the allowed region. In contrary, we cannot determine the 
value of S for a relatively short baseline like L = 295km and 1000km even if the probabilities 
are measured for two different energies, because of an almost identical overlap of regions 
with the same color for the two different energies (left and right figures). The overlapping 
is over the whole range of angles 0° < S < 360°. We can determine the slope of these two 
regions as about —0.006, which is almost equal to the coefficient of e, namely —1/170 in 
(4.6). From this observation, we conclude that the uncertainty of Am|]^ of more than 0.6% 
prevents us from determining S from f^ -^ u^ oscillations only in relatively short baselines. 
Or in other words, for the case of short baseline length the determination of the CP phase 
6 will become possible, if the uncertainty can be decreased below 0.6%, as demonstrated 
in the following section. 

5. CP Sensitivity with More Small Uncertainty 

In this section, we discuss the possibility for diminishing of parameter uncertainties includ- 
ing Arrig]^ in order to measure 6. We also investigate how the CP sensitivity is improved 
in a relatively short baseline superbeam experiment after diminishing of parameter uncer- 
tainties. It is discussed that the uncertainty of Am|]^ can be reduced up to 3% at 3-cj 
by loading 0.1% Gd in the detector of SK [39] and that the uncertainty of sin^ 612 can be 
reduced up to 2% at l-o" by using the reactor experiment with the baseline of L ~ 60km 
[40]. From these analysis, there is a possibility for diminishing the uncertainties up to 1% 
at 1-cj C.L. for both Am^i and 9i2 in near future experiments. The uncertainty of Amg^ 
can be reduced up to 1% at the T2K experiment and the NOi^A experiment. However, it 
is required that the uncertainty has to be reduced one more order of magnitude in order to 
receive the sensitivity for the CP phase as discussed in the previous section. So, we consider 
two baselines. We measure the value of ^m\i as precisely as possible in one baseline and 
then measure the value of 6 by using only v^ -^ v^ oscillations in another baseline. As 
one of the real models, we consider the T2KK experiment [19]. In this experiment, the 
neutrinos emitted from Tokai are observed in two detectors at Kamioka and Korea. One 
of the aims in the T2KK experiment is to decrease the systematic error by using the same 
beam, but we use this experiment to reduce the uncertainty of ^m\i. Here, we assume 
that we perform the precise measurement of Am^^ in Korea WC detector at L = 1000km 
with 500kt fiducial mass and then we measure the value of 5 in Kamioka WC detector at 
L = 295km with the same size. The size of detectors assumed here is larger than those 
considered in [19]. 

At first, we set the true value of Amgj^ as 2.5 • 10~^eV^ and calculate A^^ as the 
function of test value at L = 1000km. Here, we consider the case of normal hierarchy. In 
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Figure 10: Precision of ^rn'^i. Blue and pink lines correspond to the case of Stme 
Ax^ is calculated considering both systematics and parameter uncertainties. 



90° and 180° 



the case of inverted hierarchy, we obtain the similar features. In our analysis, we consider 
6 as free parameter and set the uncertainties of 1% for Atti^i and Ou- We use the same 
uncertainties as in sec. 3 for other parameters. Figure 10 shows the value of A^^- Pink and 
blue lines correspond to the case of 6true = 90° and 180°. We take into account not only 
the systematic error but also the uncertainties of various parameters in calculating Ax^- 
From this figure, the test values separated from the true value more than 0.2-0.3% can be 
excluded at l-cr C.L. with 2 d.o.f. 

Figure 11 shows the value of Ax^ as the function of test value 6, fixing the uncertainty 
of A?n|^ to 0.2% at 1-a. About uncertainties of other parameters, we use the same as in 
figure 10. 

We found that there is the CP sensitivity even at relatively short baseline like L = 
295km due to the decrease of parameter uncertainties. For 5true = 90°, the allowed range 
is 86° (190°) at 1-a {2-a) and for Stme = 180° the allowed range is 100° (148°) at 1-a {2-a). 
In such short baseline, we have an advantage in the statistical point of view. Therefore, this 
strategy may be better if we can determine the parameters precisely in future experiments. 

In figure 12, we show the ^13 dependence of the CP sensitivity in the baseline L = 
295km when 5{Am'^^) / Am'^'^ is fixed at 0.2%. The experimental setup is the same as in 
sec. 3. 

From figure 12, we can see that the CP sensitivity becomes worse mildly according to 
the decrease of ^13 as in figure 6. In the case of 5true = 90° and 180°, the allowed range is 
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Figure 11: Ax^ by assuming the true values, 6 — 90°, 180°. Pink and blue lines represent the 
Ax^ considering only systematics and including also parameter uncertainties. 
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Figure 12: 6113 dependence of CP sensitivity. In left and right figures, we take Stme — 90° and 
180°. Red, green and blue lines show 1,2 and 3-<t C.L. lines respectively. 



about 222° and 160° respectively at 1-a for the value of sin'^ 2^13 = 0.03. 

Finally, let us comment the dependence of CP sensitivity on the systematic errors. In 
the baseline of L = 295km, the first term Ai in eq.(2.11) can be neglected and only the 
second term A2 contributes to the determination of 6. Namely, we obtain the informa- 
tion of 6 by the energy dependence of QE f^ events. In our analysis, we take the signal 
normalization of QE u^ events as free from the beginning. Therefore, the signal normal- 
ization hardly affect the CP sensitivity in this baseline. This is in contrast to the case of 
L = 5000km. 
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6. Summary and Discussion 

In this paper, we have explored new possibihties for experiments to be performed in the 
next decade for the case that sin^ 2^13 is larger than 0.03. The ^13 will be determined in 
the next generation reactor experiments. Concretely, we have investigated whether the CP 
phase can be measured by only f^ -^ f^ oscillations independently of v^ -^ Ve oscillations. 
If we can measure the CP phase in two different channels independently and there is a 
difference for these values, this would be considered as an evidence of new physics. Below, 
the results obtained in this paper are listed. 

• At first, we have investigated the energy and the baseline where the cos 5 term in- 
cluded in Pfj_^ becomes large by using both numerical and analytical methods. As the 
result, we found from (2.17) and (2.18) that the coefficient A^^ has its largest value 
around E = 0.43GeV and L = 5000 and 10000km in the earth mantle. The difference 
of the probabilities attains the maximal value about 0.2 due to the CP phase effect 
even after the averaging over the energy. 

• Next, we have considered the same beam and the detector as in the T2HK experiment 
but the baseline L = 295km and 5000km and have calculated the CP sensitivity by 
using x^ method. As the result, when the uncertainty of Am|^ is larger than 1%, the 
CP sensitivity has its best value for the baseline around L = 5000km. The allowed 
range becomes 50° (88°) at l-a in L = 5000km when Stme = 90° (180°). On the 
other hand, we have almost no sensitivity for the CP phase in the case of L = 295km 
because of the effect of parameter uncertainties. 

• We have shown both numerically and analytically that the uncertainty of Arrigj^ is 
particularly important in determining the value of 5. If the uncertainty of Am|]^ will 
be less than 0.6%, a certain sensitivity for 6 will be obtained even in the relatively 
short baseline like L < 1000km. 

• We have explored the possibility for measuring the CP phase in relatively short 
baseline by diminishing the uncertainty of Arrig^. As concrete experimental setup, 
we have considered the T2KK experiment. We used one baseline from Tokai to 
Korea in order to determine A?ti|]^ precisely and used the other baseline from Tokai 
to Kamioka to measure the CP phase by t'^ -^ t'^ oscillations. As the result, the 
uncertainty of Am|]^ can be reduced up to about 0.2% and the allowed range becomes 
about 222° (160°) at 1-cr when 6true = 90° (180°) at 1-a. So, we have the possibility 
for measuring the CP phase although the sensitivity is not so good. 

In future, the mixing angles and the mass squared differences are precisely measured 
in various kinds of experiments. If ^13 is found in the next generation reactor experiments 
in addition to this improvement, it will be very important to consider the strategy for 
exploring the new physics. As one of the strategies, the possibility for measuring the CP 
phase in v^ — > f^ oscillation is interesting. We will further investigate in our next work 
how the contribution of new physics appears in future experiments. 
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A. Derivation of Approximate Formula for A^^ 

In this appendix, let us derive the approximate formula of A^^ used in sec. 2 for the case 
of constant matter. 

At first, in the case of ^23 = 45° and the symmetric matter profile, the probability is 
given by 

^MM = \^^ cos 6 + C^^ + D cos 26, (A. 1) 

from the exact formula [30, 32]. Here, we can neglect D because of the higher order terms 
of small parameters. Concrete expression for A^^ [32] is given by 

cvclic 

A,, = Y. 8[J.A.,A3,A..(A.-A30-(^..).] ^^^ ^,^. ^.^ ^, ^ ^.^ ^, ^^ (^2) 

(ijk) ^jk^ki 



where 



{A^^)k = Ail Jr X [A3iAfe(c?2 - s?2) + Ais?2 " ^IA2] ■ (A.3) 



In the energy region around E = IGeV, we have Ai = 0(A2i),A2 = 0(A2i) and 
A3 ~ A31. From these order counting, we can neglect Ai and A2 compared with A3. On 
the other hand, we cannot neglect A2X included in the oscillating term for long baseline 
such that the relation A21 = Am2iL/4:E = 0(1) is satisfied. If we leave only leading order 
terms in small parameters A21 and S13, the expression for A^^ is reduced to 

, 8J^A2iA3i[Ai(Ai - A31) + A3iA2icf2] x/ • X' ■ a' 
^MM - X2~7^ cosA23smA3ismAi2 

^31^21 



H ~„ -„ cosA3ismAi2smA23 



8J^A2iA3i[A2(A2 - A31) + A31A21CI2] 
4J^A2iA3i[Ai(Ai - A31) + A31A21C22] 



A2 A2 



sin A'lalsin A21 + sin(A3i - A23)} 



^ 4J.A2lA3l[A2(A2 - A31) + A3lA2lcf2] . A' / • A' ■ tA' A' M 

+ T2 a2 smAi2{smA2i - sm(A3i - A23)} 

^32 

sin An 



A^ A^ 

^32^21 



4J^A2lA3i . ~,2 
sin r- 

A2 



21 

21 

Ai(Ai - A31) + A3iA2lC?2 ^ A2(A2 - A31) + A31A21C22 



A2 A2 

^31 ^32 



4J^A2iA3i 

A2 



^Jr^21^31 . X/ ■ /oA' A' ^ 

H -^ sm A12 sin(2A3i - A21) 
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Ai(Ai - A31) + A3iA2lcf2 A2(A2 - A31) + A3iA2lC^2 
A2 A2 

^i^ sin A'l, [Ai - 2A2ic?2 + A2] - ^i^ sin A'^i sin(2A;,i - A',,) 
A21 A21 

^ — sm A21 [a - A21 cos 26*12] r 

A^i A21 



sin Afi [a - A21 cos 26'i2] ^^ sin A21 sin(2A32). (A.4) 



This is the derivation of (2.11). In this approximate formula, the MSW effect due to the 
1-3 mixing angle is not included, so there is some differences from the exact one in high 
energy. On the other hand, this approximate formula coincides well in low energy. In this 
expression, if L/E is small enough, the first term can be neglected and is reduced to 

A^^ ~ _li^sinA'2isin(2A;52) - -^i^ sin A'21 sin(A;3i) cos(A;3i). (A.5) 
A21 A21 

This is the third term of (4.1). 
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